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A feasibility study is presented to precisely measure the masses of charged sleptons and the 
lightest neutralino in e+e~ f^Rfj-R and e+e~ fifi production. Different mSUGRA 
^ '. scenarios with small slepton-neutralino mass differences of a few GeV are invesigated. The 

CN I analysis is based on a detailed simulation of the final state lepton energy spectra including 

background reactions and assuming realistic experimental conditions at the Tesla e"'"e~ 
, Linear Collider. Effects on the mass resolutions due to beams colliding head-on or under a 

O I crossing angle are discussed. 

&: 1 Introduction 
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In the minimal supersymmetric standard model the lightest superparticle (LSP), which is stable 
in _R-parity conserving scenarios, is a natural candidate for cold dark matter. Measurements 
5^ ! of the Wmap experiment [1] have determined the amount of cold dark matter in the universe 
to be VtcDKih'^ = 0.113 ± 0.009. As a consequence of this result, considerable constraints 
on the parameters of supersymmetric theories can be set [2]. For example, in minimal super- 
gravity (mSUGRA), where usually the LSP is the neutralino Xi, the favoured particle spectra 
are typically characterised by small mass differences between the lightest slepton (usually the 
stau) and neutralino (co annihilation region) or between the light chargino and neutralino (focus 
point region). Collider experiments will be essential to test the supersymmetry dark matter hy- 
pothesis. But only a future e^e~ Linear Collider will be able to provide precise enough mass 
measurements of the light SuSY particles in order to compute the relic density of the LSP to an 
accuracy comparable with Wmap and other proposed experiments. 

Experimentally challenging is the coannihilation region where the slepton-neutralino mass 
difference Am may be as small as a few GeV. In this note we investigate the capabilities of a 
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TeV e+e" Linear Collider like Tesla [3] to detect and measure jlnjlR and fifi pair production 
in models with small Am in order to determine precisely the masses of the light sleptons as well 
as of the neutralino x'( from the energy spectra of the final states. The case studies are based 
on spectra from a SPS la inspired model [4], where the neutralino mass has been increased to 
accommodate a small mass difference Am and from the benchmark model D' proposed in [5]. 

The most severe background comes from four fermion production via two-photon (77) pro- 
cesses, e+e~ e~^e~ ff, and it is important to have an efficient veto of the spectator elec- 
trons/positrons down to very small scattering angles. The present study also addresses the con- 
sequences for background suppression in head-on collisions, as envisaged in the TESLA design, 
versus the option of having a crossing angle of 2 ■ 10 mrad between the colliding beams. 



2 Event generation 

The processes under study are the pair production of scalar muons and scalar taus 

eL^fl fj'RfJ'R At"^X?Ai"X?, (1) 
et^R ^ hh ^ r+Xir~Xi- (2) 

The beams are assumed to be longitudinally polarised, right-handed electrons and left- 
handed positrons e^, which considerably increase the production cross sections and suppress 
background. 

The slepton decay i~ iiito ^^i ordinary lepton and neutralino is isotropic and pro- 

duces a flat lepton energy spectrum. The endpoint energies £'+ /_ are related to the primary 
slepton and the neutralino masses (neglectging the lepton mass) 



m- - .rs^^^ (4) 

m, = (5) 

One observes that with decreasing slepton-neutralino mass difference Am the endpoint energies 
get lower and come closer together. For r decays the spectrum is quite distorted, the maximal 
energy of the observable particles coincides with the upper edge but the lower edge at 
is completely diluted. However, the shape of the energy distribution can still be used to extract 
the f mass with high precision. 

Events are generated with the program Pythia 6.2 [6] which includes beam polarisations 
Pe±> initial and final state QED radiation as well as beamstrahlung a la CiRCE [7]. The decays 
of r leptons are treated by Tauola [8]. 
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Figure 1: Energy deposits [GeV/cm^] from beamstrahlung at ^/s = 500 GeV for head-on 
collisions (left) and crossed beams (right) in a plane perpendicular to the beams at 3.7 m from 
the interaction point [10] 

The detector simulation is based on the detector proposed in the Tesla TDR [3] and imple- 
mented in the Monte Carlo program SiMDET 4.02 [9]. The main detector features are excel- 
lent particle identification and measurement for a polar angle acceptance 9 > 125 mrad. The 
forward region^ is equipped with electromagnetic calorimeters in order to effectively veto elec- 
trons and photons down to very small scattering angles. For 9 < 27.5 mrad a highly segmented 
LCAL calorimeter situated at a distance of 3.7 m from the interaction point has recently been 
proposed [10]. For head-on collisions, as envisaged at TESLA, a beam pipe of radius r = 1.2 cm 
appears feasible and the calorimetric coverage may start at 3.5 mrad (assuming 1 mm clear- 
ance). For beams colliding under a crossing angle of 2 • 10 mrad there are two dead regions 
and the exiting, disrupted beams require a larger beampipe of r = 2.0 cm, giving a minimal 
tagging angle of 5.7 mrad. In both cases e, 7 tagging at low angles is very demanding above 
the huge background of energy deposits caused by beamstrahlung, as shown in figure 2 [10]. 
The beamstrahlung halo varies rapidly around the beam pipes. Further, the angle between the 
magnetic field and the crossed beams leads to an enhanced (factor 2) and asymmetric energy 
distribution. An algorithm based on topological properties of electromagnetic showers is used 
to veto electrons and photons a few mm close to the beam pipe for energies above 50 GeV with 
high efficiency [11] in both designs. 

For Tesla a geometry with crossed beams is an option which may facilitate beam diagnos- 
tics. A normal conducting linear collider like Nlc/Glc has to be operated with a crossing angle 
and the short bunch separation may lead to additional, disturbing pile up effects, depending on 
the calorimeter read-out. This difficulty will not be further discussed. 

'The acceptance at an angle 9 in the forward region is symmetric in the backward direction at vr — 
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3 Case study SPS la inspired model 



The aim of this paper is to investigate the experimental problems related to small lepton- 
neutralino mass differences and to compare the results with more conventional benchmark spec- 
tra. As reference serves the mSUGRA scenario SPS la where detailed studies of slepton pro- 
duction in the continuum [12] and at threshold [13] exist. The case study is based on a SPS la 
inspired model, where the masses of the the light sleptons are kept unchanged, m^^ = 143 GeV 
and nif-^ = 133 GeV, and the neutralino masses are increased. Within mSUGRA a fi — mass 
splitting Anif = 8 GeV can e.g. be accomplished by shifting the common gaugino mass mg 
to 70 GeV and the common scalar mass ??7i/2 to 320 GeV and leaving the other parameters as 
tan/5 = 10, Aq = -100 GeV and sign fi>0. 

The strategy will be to determine the neutralino mass from the lepton spectra of jlRjlR and/or 
eRS-R production (eqs. (5) and (4)) and then use this value as input in the r decay spectra of the 
f analysis. Optimal collider operating conditions are chosen such as to maximise the signal 
cross section and to minimise the background: a center of mass energy of ^/s = 400 GeV 
and beam polarisations of Ve- = +0.8 for right-handed electrons and of Ve+ = —0.6 for left- 
handed positrons. Other SuSY processes are kinematically not accessible (below production 
thresholds) or strongly suppressed (fif2, crCl and XiX^)- The SM contribution from WW 
production is negligible due to the beam polarisation and efficient selection criteria. The es- 
sential remaining background comes from two-photon processes e+e" — *■ e~^e~i~^i~ with cross 
sections of (9(10® fb). The simulations assume an integrated luminosity of C = 200 fb~^, 
corresponding to about one year of data taking. 

3.1 Detection of e+e~ — ^ i^ni^R 

For the purpose of illustration a small smuon-neutralino mass difference of Am = 8 GeV is 
assumed. Such a situation may occur for low values of tan (3 with the light smuon and stau 
being almost degenerate. The masses used in the simulation of e^e^ — > /i^/i^ /^"""X? 
are m^^ = 143.0 GeV and m^o = 135.0 GeV. The relevant signal and background cross 
sections are 0-/!^^^ = 120 fb, aww = 1000 fb and Uee^^ = 9.4 ■ 10^ fb. 

The following event selection criteria are applied: (1) veto of any electromagnetic energy 
(e or 7) above 5 GeV in the forward region 6 < 125 mrad, (2) two muons within the polar 
angle acceptance —0.90 < Q^cosO^ < 0.75, (3) acoplanarity angle A^'^^ < 160°, (4) missing 
momentum vector inside active detector cos6'p-^.^^ < 0.9, (5) muon energy > 2 GeV, (6) 
transverse momentum of di-muon system p'^ > 5 GeV. The overall signal efficiency is 60 %, 
while the remaining background is negligible. The most effective cuts against WW production 
are (2) and (3). The 77 /i/x background is reduced by a factor of ~ 10~®, particularly 
powerful are cuts (3) and (6). 

The resulting spectra of the muon energy i?^ are shown in figure 2. The rectangular shape 
with the steeply rising edges is clearly observable. The background is roughly a factor of 2 
larger in a configuration with beam crossing angle compared with head-on collisions. However, 
the influence on the endpoint energy measurements are in both cases negligible. A fit to the 
spectrum yields the endpoint energies E_ = 3.270±0.015 GeVandE+ = 18.480±0.39 GeV, 
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Figure 2: Muon energy spectra from the reactions e^e^ij — > I^rI^r — > A^^Xi Xi (points) 
and e^e^ e+e~ yU^/x^ (histogram) without (left) and with (right) beam crossing angle. 
SPS la inspired model, m^^ = 143 GeV, Am = 8 GeV, ^ = 400 GeV and £ = 200 fb"^ 



from which the strongly correlated masses of fiR and Xi are derived with an accuracy of 5m^^ = 
0.18 GeV and 5m^o = 0.17 GeV, see table 1. 

The precision on the neutralino mass can be considerably improved by analysing the se- 
lectron production e^e^ — e^e^ — e^Xi ^~Xv The cross section is roughly a factor of four 
higher while the two-photon background remains low. From the results of table 1 one concludes 
that a neutralino mass resolution of 6m^a = 0.08 GeV is achievable, which can be further im- 
proved by a combined analysis of both slepton production processes. 



e+e- ^ jlnjlR m^^ = 143.0 ± 0.18 GeV m^o = 135.0 ± 0.17 GeV 
e+e" ^ crCr mg^ = 143.0 ± 0.09 GeV m^o = 135.0 ± 0.08 GeV 

Table 1: Expected mass resolutions for jlR, cr and Xi from lepton energy spectra in a SPS la 
inspired scenario assuming C = 200 fb~^ at ^/s = 400 GeV 

Compared with the original SPS la scenario [12], the mass resolutions obtained from the 
first and second generation slepton pair production do not degrade with decreasing slepton- 
neutralino mass difference down to values of Am = 8 GeV, and possibly even lower to 5 GeV. 



3.2 Detection of e+e f\f\ 

For the simulation of e^e^ ^t^i a stau-neutralino mass difference of Am = 

8 GeV is assumed, leading to m^^ = 133.2 GeV and m^o = 125.2 GeV. The cross section 
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is (Tf^fj = 140 fb, SM background considered have cross sections of aww = 1000 fb and 
(TeerT = 4.5 ■ 10^ fb. For the identification and reconstruction of r pairs, one tau is required to 
decay into hadrons while the other may decay hadronically or leptonically. The measurement 
of energy spectra is restricted to the hadronic decay modes r (11.1%), r pUr 

TT^TT^z/,- (25.4%) and r SnVr n^n^n^Vr + h^h^h^Ut (19.4%). A difficuky may pose 
the r polarisation from fi ry^ decays, which has an influence on the pion energy spectrum, 
but does not affect the shape of the p and Svr spectra [14]. The r polarisation can either be taken 
from calculations or can be determined simultaneously [12], if one wants to keep the pion decay 
channel. The leptonic 3-body decays r eu^Ur^ t P^f^^r exhibit a very weak dependence 
on uif and are therefore discarded in the analysis. 

The following criteria are applied to select fifi events: (1) veto of any electromagnetic 
energy (e or 7) above 5 GeV in the forward region 6^ < 125 mrad, (2) two r candidates within 
the polar angle —0.75 < cos6't- < 0.75, (3) invariant mass of r jet m,- < 2 GeV, (4) tau 
energy 2 < i?,- < 25 GeV, (5) acoplanarity angle A^^"^ < 160°, (6) missing momentum 
vector cos6'p>^-^^ < 0.8, (7) transverse momentum of di-tau system 3 < pY < 25 GeV, (8) 
combined cut on ^ and /S.cjf'^ . An effective reduction of 77 background is possible by 

using observables defined in the plane perpendicular to the beams. Cuts (5) and (7) are not as 
efficient as in the smuon analysis. They have to be supplemented by criterion (8), illustrated in 
figure 3 as correlations between ^ ^, the sum of r momenta projected onto the transverse 

thrust axis Ti, and the acoplanarity angle A^"""". A cut as indicated by the curves has little 
influence on the signal and reduces the background by factors of order 50. The residual fifi 
efficiency is ~ 18%. The two-photon background is suppressed by ~ 1.5 ■ 10"^ for head-on 
collisions, respectively ~ 2.5 ■ 10^^ in the case of crossed beams. Contributions from WW 
production are negligible. 




Figure 3: Correlations of ^ ^ wrt thrust axis T±_ versus acoplanarity angle ^(jf^ in a plane 

perpendicular to the beams for e^e^ ^ n fi ^ t^Xi t^Xi (left) and e^e^ e~^e~T^ 
(right). The curves indicate the separation between SuSY and 77 processes. SPS la inspired 
model, mf, = 133 GeV, Am = 8 GeV, = 400 GeV 
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The energy spectra of i?^, Ep and E'^^, are shown in figure 4 for collisions with crossed 
beams. Clean signals above a very small and controllable background are observed. The spectra 
are fitted to full simulations with different fi mass hypothesis, yielding a sensitivity of ^m^^ = 
0.14 GeV © Sm^o for the Et^ spectrum, and 5mf-^ = 0.10 GeV © Sm^o for the Ep and E-s-,, 
spectra, each. Assuming a conservative estimate of Sm^^o = 0.10 GeV one expects a combined 
uncertainty of 5mf-^ =0.14 GeV for the the stau mass. A lower background expected in head-on 
collisions or a twofold higher rate (accounting for additional two-photon processes not covered) 
do not change the resolution significantly. 
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Figure 4: Hadron energy spectra E.,, of r ^ nVr, Ep of r ^ pu^ and E^^t, of r ^ Svrz/^ decays 
from the reactions e^e^ fi ^ t^Xi ^^id e\e~^ — > e^e'r^r^ assuming crossed 

beams collision. SPS la inspired model, nif^ = 133.2 GeV, Am = 8 GeV, ^/s = 400 GeV 
and C = 200 fb~^ 

It is interesting to note that the precision achieved in the present study is about twice as 
good as for the original SPS la scenario [12]. The reason is that for small stau-neutralino mass 
differences the energy spectra are narrower and are falling steeper at the edges. In particular, 
a determination of the high endpoint energy, eq. (3), is therefore more precise, which directly 
enters the mass calculation. 

Examples of a further reduction of the stau-neutralino mass difference to Am = 5 GeV and 
Am = 3 GeV are shown for the p energy spectra in figure 5. The overall event rates decrease 
but are still sufficient to clearly determine the high endpoint energy. Including mass effects, the 
maximal r energy and its relation to the fi mass are given by 

2 2 

m| — mlo 

^+ = ^ ^ ^ith /. = 1 - 0.5 (m^./Am)^ (6) 

6mf^ ~ 0.38 SE+ © 6m^o . (7) 

A very simple fit of a straight line to the upper part of the spectra yields a maximal energy 
of E+ = 12.94 ± 0.52 GeV which translates to 6mf, = 0.22 GeV for a mass difference of 
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Figure 5: Hadron energy spectra Ep of t pUr from eje^ t^iTi ^ t'^Xi t Xi with 
Am = 5 GeV (left) and Am = 3 GeV (right) assuming head-on collision. SPS la inspired 
model, ruf^ = 133.2 GeV, = 400 GeV and C = 200 fb~^ 

Am = 5 GeV. For the scenario with Am = 3 GeV one obtains = 7.60 ± 0.72 GeV or 
5mf-^ = 0.28 GeV. These resolutions for very small stau-neutralino mass differences are very 
encouraging and can be further improved with more refined spectra analyses. All results for the 
fi mass precision obtained under different assumptions are summarised in table 2. 



4 Case study model D' 

The post-WMAP model D' [5] (parameters mo = 101 GeV, mi/2 = 525 GeV, tan (3 = 10, 
Aq = GeV, sigHyU < 0) is a typical benchmark in the coannihilation region. The relevant 
masses are the light sleptons m^^ = mg^ = 223.9 GeV and m^^ = 217.5 GeV and the 
neutralino m^o = 212.4 GeV. For the analysis the collider is chosen to operate at ^/s = 
600 GeV with polarised beams of V^- = +0.8 and Ve+ = —0.6, which essentially suppresses 
all other Susy production processes. The luminosity is assumed to be £ = 300 fb^\ 

From the discussion of section 3.1 it is safe to conclude that the neutralino mass can be 
determined to an accuracy of Sm^o ~ 0.1 GeV, which will be assumed in the following stau 
analysis. 



4.1 Analysis of e+e T{'f^ 

The study of e^e^ 'T^iT^i ~^ '^'^Xi t~Xi proceeds along the same lines as described in 
section 3.2. The cross section is cTf^f^ = 50 fb. The event selection is identical except a tighter 
cut (6) on the tau energy of 2 < < 15 GeV, due to the small stau-neutralino mass difference 
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Am = 5.1 GeV. The resulting overall efficiency is 7.6%. In addition to e'^e" e'^e''TT 
a hadronic background from two-photon quark production 77 qq is included. The major 
contribution comes from charm quarks where the final states mimic a tau pair. There is still 
an uncertainty on some missing graphs concerning real or 'direct' 77 interactions, which may 
contribute with a similar magnitude. Fortunately, the total background can be kept small at a few 
percent level, and an increase by a factor of ~ 1.5 would not change the results significantly. 
In a final experiment the 77 background will be precisely measurable and should be reliably 
calculable. 
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Figure 6: Hadron energy spectra of t ^ ixu^-, Ep of t ^ pu^- and E^^, of r — Stti^t- decays 
from the reaction e^e^ — ^ fi — > t~^Xi two-photon production assuming head-on 

coUision. Model D', m^^ = 217.5 GeV, Am = 5.1 GeV, = 600 GeV and C = 300 fb"^ 



The energy spectra of E^,, Ep and E^^, are shown in figure 6 for head-on collisions. One 
observes narrow distributions with a small background increasing towards low energies. From 
fits to the spectra assuming various fi mass hypothesis, one obtains uncertainties in the mass 
determination of 5mf^ = 0.19 GeV © Sm^o for the Et, spectrum, Snir-^ =0.12 GeV © Sm^o 
for the Ep spectrum and dm^j = 0.14 GeV © Sm^o for the E-^^, spectrum. Combining the 
information of all measurements and folding in the neutralino mass error {Sm^o = 0.1 GeV) 
the stau mass can be determined as m^^ = 217.5±0.15 GeV. Despite the lower event rates the 
achievable accuracy is very similar to that obtained for the SPS la inspired model (see section 
3.2). This supports the observation that narrow distributions have a somewhat higher analysing 
power. The results of the simulation of model D' are given in table 2. 



5 Conclusions 



Simulations of slepton production e~^e PrP-r and e"*"e rifi in scenarios with small 
slepton-neutralino mass differences Am are presented under realistic experimental conditions. 
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scenano 



mrAGeY] ^m^, [GeV] Am [GeV] dQcDMh'^ 



model D' 



SPS la inspired 



133.2 0.14 8 1.7% 

0.22 5 4.1% 

0.28 3 6.7% 

217.5 0.15 5.1 1.9% 



Table 2: Expected accuracies on the stau mass 6mf-^ and the resulting precision on dark matter 
densitiy 5 VtcDNih'^ for variants of a SPS la inspired scenario {C = 200 fb"^ @ 400 GeV) and 
the benchmark model D' (C = 300 fb"^ @ 600 GeV) 

Choosing the proper collider energy and beam polarisations, the decay lepton energy spectra 
are analysed and allow the masses of the light sleptons and the neutralino Xi to be determined 
very precisely, to an accuracy of one per mil. The results of various case studies with ranges 
Am = 8 GeV 3 GeV are summarised in table 1 and table 2. 

Methods are devised to efficiently suppress the most serious background from 77 four- 
fermion production e^e~ e'^e^ ff, with f = r, q, which becomes more important with 
decreasing Am. It is mandatory to have excellent veto capabilities for the scattered e"^ down 
to very small angles as close as possible to the beam pipe. Low angle tagging can be more 
easily achieved in head-on collisions of the Tesla design. In the case of a crossing angle of 
2 ■ 10 mrad the 77 background is larger by a factor of 1.5 — 2. However, one does not expect 
significant degradations of the mass resolutions as long as the background can be kept small. 



For the model D', which is typical for the coannihilation region, the achievable f mass 
resolution from energy spectra is (5mf^ = 0.15 GeV. The present analysis can be compared 
with an alternative method, based on a cross section measurement of fifi production close to 
threshold [15]. Under favourable conditions and assuming a luminosity of 500 fb^^ an accuracy 
of 0.54 GeV is quoted, which is considerably worse than the value obtained from energy spectra 
of the decay particles. 

It is interesting to put the results in the context of a dark matter scenario. Within the coan- 
nihilation region the dark matter content of the universe is essentially determined by Am, re- 
spectively the masses rrif^ of the stau and m^o of the neutralino. A calculation of the relic dark 
matter density with the program MICROMEGAs [16] and using the uncertainties of the present 
study shows that VtcDMh"^ can be predicted at the few percent level, see table 2. This emphasises 
the important role of a future TeV Linear Collider to provide high precision measurements. 
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say < 0(10%). 
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